Spatial and Spectral Multifractality of the Local Density of States at the Mobility 

Edge 
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We performed numerical calculations of the local density of states (LDOS) at disorder induced 
localization-delocalization transitions. The LDOS defines a spatial measure for fixed energy and a 
spectral measure for fixed position. At the mobility edge both measures are multifractal and their 
generalized dimensions D(q) and D(q) are found to be proportional: D(q) = dD(q), where d is the 
dimension of the system. This observation is consistent with the identification of the frequency- 
dependent length scale L u cx oj _1//ci as an effective system size. The calculations are performed for 
two- and three-dimensional dynamical network models with local time evolution operators. The 
energy dependence of the LDOS is obtained from the time evolution of the local wavefunction 
amplitude of a wave packet, providing a numerically efficient way to obtain information about the 
multifractal exponents of the system. 
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In the presence of disorder electronic states in phase- 
coherent systems can get localized. For dimensions d > 2 
such a disordered system can exhibit a transition from lo- 
calized to extended eigenstates as a function of energy or 
strength of the disorder, similar to continuous thermo- 
dynamic phase transitions H. This transition is accom- 
panied by a diverging length scale £(E) oc \E — E c \~ y 
that can be identified for localized states with the lo- 
calization length. Here E c is the critical energy of the 
transition. On length scales shorter than £ (E) the eigen- 
states exhibit strong fluctuations Bj . Sufficiently close to 
the transition the localization length exceeds the system 
size so that the wave functions fluctuate on all length 
scales up to the system size L. The squared modulus 
of these critical eigenstates forms a multifractal measure 
p(r) = \ip(r)\ 2 Bj. The scaling behavior of multifractal 
measures is determined by an infinite set of exponents 
D(q), called generalized dimensions. The multifractal 
fluctuations lead to anomalous behavior of the diffusion 
coefficient f|,||]. 

The localization properties of the wavefunction are also 
reflected in the spectral properties of the system. Local- 
ized states correspond to a pure point like spectrum, ex- 
tended states to an absolutely continuous spectrum, and 
critical states to a singular continuous spectrum. The 
latter kind of spectrum was observed in quasiperiodic 
systems Q, systems on incommensurate structures [Q, 
in crystals in a magnetic field and also in random 
one-dimensional systems || . For these systems the spec- 
tral measure excited by a wave packet is multifractal. Of 
particular interest has been the implications of the lo- 
cal spectra for dynamical properties [ p~0jJ5~l| ] . Ketzmerick, 
Petschel, and Geisel ]l(J found that the return probabil- 
ity of a wave packet decays in time t asymptotically as 



t D ( 2 \ where D(2) is the correlation dimension of the 
associated spectral measure, the local density of states 
(LDOS). 

The two aspects of spatial and spectral multifractality 
were brought together by Huckestein and Schweitzer who 
showed that enhanced return probability of wavepackets 
at the mobility edge of quantum Hall systems could be in- 
terpreted both from the spectral Gq] as well as the spatial 
properties of the local density of states || . They showed 
that the generalized dimensions D(2) and D(2) charac- 
terizing the second moments of the spatial and spectral 
measures, respectively, are related by the spatial dimen- 
sion d = 2 of the system, D(2) = dD(2). This follows 
from dynamical scaling by which the two-particle spec- 
tral function S(q, u>) becomes a function of qL^ at crit- 
ically Q. Here, L u = (p(E c )uu)~ 1 / d can be interpreted 
as the system size with mean level spacing u> and p(E c ) 
is the density of states at the mobility edge The 
length scale introduced by the finite frequency w cuts 
off correlations and acts as an effective system size [[l4} . 
Recently, the same relation was shown to hold in three 
dimensional systems at the Anderson transition fHI . 

In this letter, we show that at the mobility edge the 
spatial and spectral structures of the LDOS are inti- 
mately related: their respective generalized dimensions 
D{q) and D(q) are proportional, 



D(q) = dD(q), 



(1) 



where d is the space dimension of the system. This re- 
lation, generalizing the result of ref. || to arbitrary q, is 
a consequence of the length scale acting quite gener- 
ally as an effective system size, not only for correlation 
functions but also for moments of the LDOS. We test the 
relation numerically for two- and three-dimensional net- 
work models |16|,|l7||. Introducing a time evolution into 
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these models [Q allows us to calculate the LDOS by 
studying the time evolution of wave packets. This pro- 
vides a new efficient method to obtain the multifractal 
exponents of the LDOS without the need to diagonalize 
large matrices. 

Before motivating the relation (|l]) let us define our mul- 
tifractal analysis. We study a normalized multifractal 
density p(x) on a d-dimcnsional domain Q. Here, Q is the 
two- or three-dimensional space or the one-dimensional 
energy axis, depending on whether we study the spatial 
or spectral aspect of the LDOS. We study the scaling of 
the box probabilities 

Pi(l) = I d d xp{x) (2) 

with respect to the box sizes /. f2j(/) is the volume of 
the ith of Ni = (L/l) d non-overlapping hypercubes of 
linear size I covering the whole domain f2 of linear size 
L. For a multifractal density the averaged gth moments, 
(P q (l,L)) = Nf 1 EiSi (-Pi(O) 9 . show power law depen- 
dence on l/L <C 1 for all real q, 

<P«(Z,£))oc(^J . (3) 

The distinguishing feature of a multifractal is that the 
exponents r(q) = (q— l)D(q) are a non-linear function of 
q. A multifractal is thus described by a non-countable set 
of exponents D(q). If the multifractal measure is taken 
from a statistical ensemble, the system average over the 
box probabilities can be replaced by an ensemble average 
for a single box. 

The LDOS is given by p(r, £)=E a S(E-E a )\<f> a (r)\ 2 , 
where the sum runs over the quantum numbers a with 
corresponding eigenstates <f> a . At energies E close to the 
mobility edge E c the LDOS exhibits multifractal proper- 
ties on length scales between microscopic lengths, like the 
lattice constant a or the elastic mean free path, and the 
system size L or the correlation length £(E). The cor- 
responding lower and upper energy scales are the mean 
level spacing A = p(E c )^ 1 L^ d and Eg oc L~ x l v , respec- 
tively. For energies u> less than Eg the correlation length 
exceeds the system size L. 

We now argue that the existence of a single frequency- 
dependent length scale implies the relation (|l|) . Con- 
sider box probabilities of the LDOS with respect to both 
space and energy 

Pi(l,uj)= [ d d x [ ° 2 dEp(r,E), (4) 
Jni(i) Je c -| 

where we consider an energy interval in the critical region 
described above. For energies u> less than the mean level 
spacing A these box probabilities scale like single eigen- 
functions, {P q {l,uj)) oc {l/L) d+T ^\ For energies larger 
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FIG. 1. Generalized dimensions D(q)/2 (O) and D(q) (□) 
for a two-dimensional network model at the quantum Hall 
critical point. Systems of size 150 x 150 (O) and 200 x 200 
(□) were used and the statistical uncertainties are smaller 
than the symbol sizes. 

than the mean level spacing, we expect the disorder av- 
erage of the moments of these box probabilities to scale 
with the same exponents but with the system size L re- 
placed by a frequency-dependent effective system size L u , 

(if \j- J • (5) 

In principle, the frequency-dependent length scale intro- 
duced in eq. (||) could depend on q. However, if we as- 
sume that there exists only one such length scale, then 
we can identify L u with (p(£' c )(jj) _1 / d by considering the 
first moment 

{Pifrw)) = p{E c )u>l° 'oc (J-^j . (6) 

The scaling of the spectral measure defined by the LDOS 
follows from eq. (^) by choosing the spatial box size I of 
the order of the microscopic length scale £, below which 
the wavefunctions become smooth functions of coordi- 
nate. The scaling of the spectral measure with respect 
to energy defines the exponents f(q) = (q — l)D(q), 

(P?(£,u)) oc {p(E c )£ d uj) 1+T(q)/d oc c 1+ ^>, (7) 

leading to eq. (|l|). 

For our computations, we consider two normalized 
measures, the spatial measure p(r) = |^ Q (r)| 2 and the 
spectral measure p(E) = A~ x p (i"o, E). The spatial mea- 
sure is normalized since the eigenfunctions are normal- 
ized, and the spectral measure is explicitly normalized 
on an energy interval E, satisfying E <C Eg. To support 
the validity of eq. (Q) we now present results of numer- 
ical calculations. Fig. (|l|) shows the functions D(q) and 
2D(q) calculated for a two-dimensional network model 
at the quantum Hall critical point |l6],[l8| . The functions 
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FIG. 2. Generalized dimensions D(q)/3 (O) and D q (□) for 
a three-dimensional network model at the Anderson transition 
in a magnetic field. Systems sizes of 20 3 (O) and 35 3 (□) wore 
used and the statistical uncertainties in the D(q) are less than 
those in the D(q)/3. 



agree within the error bars. In order to check the depen- 
dence on the spatial dimension d of the system we study 
a three-dimensional network introduced recently |l7| . In 
contrast to the two-dimensional network here a band of 
extended states appears. At the mobility edge of this sys- 
tem we again find good agreement with eq. (|l|) as seen 
in Fig. ®. 

We now describe our new numerical method of ob- 
taining the LDOS. We use two- and three-dimensional 
network models, describing the integer quantum Hall ef- 
fect [[l6| and the so called quantum Hall- insulator |L7|, 
respectively. Both models are extended to dynamical 
network models by providing them with a unitary time 
evolution operator U for discrete microscopic time steps 
[pTstl . Within these models the evolution — involving full 
quantum interference — of arbitrary initial states can be 
obtained easily by iterative application of U. In par- 
ticular, the time evolution of a state tp initially sharply 
peaked at coordinate ro yields the temporal Green's func- 
tion G(tq, Tq; t). The LDOS is then the Fourier transform 
of this quantity. This provides an efficient method of cal- 
culating the LDOS, without the need to diagonalize the 
operator U or the associated Hamiltonian. 

Starting point of our calculations is the 2D-network 
model introduced by Chalker and Coddington jl(| to 
describe non-interacting spinless electrons in the inte- 
ger Quantum Hall regime. It exhibits a localization- 
delocalization-transition characteristic for the QHE: ex- 
cept for the critical energy E c all states are localized, 
with a localization length £ cx \E — E c \~ u , v — 2.3. Ori- 
gin of the model and exact definitions can be found in 
refs. |L6| and |18[ ], here we give only a brief description 
restricted to our purposes. 

The network consists of 2 x 2 scattering matrices as 
nodes which are arranged on a square lattice and con- 
nected by one-dimensional unidirectional channels, called 



links. The scattering matrices contain coefficients ti m 
describing transitions from electron states on incoming 
links tpi to outgoing link states ip m . At the critical 
point the transmission amplitudes are of constant value 
T m i = \t m i\ 2 = 1/2, while the disorder is given by 
randomly distributed arguments of the coefficients t m j. 
States (or wave functions) ip on the network are N- 
dimensional complex vectors ip = {ipi}i=i,2,...,N, where 
the Z-th component denotes the complex amplitude on 
link I. The network operator U is defined by its action 
on single link states ei = {Sik}k=i,...,N, 



(8) 



where t m i and t n i are the transmission coefficients from 
an incoming link I into two outgoing links m and n at a 
node [pf. 

Here we use U as the time evolution operator in dis- 
crete time steps for states on the network. To motivate 
this, consider a particle at critical energy E c in the in- 
coming link I of node j, that is described by the state 
= e/. After a characteristic time r (we choose r = 1 
in the following) the incident wavepacket has passed the 
scatterer and thereby split into two outgoing packets in 
the channels m and n (see fig. 1 in ref. |l^|). This process 
corresponds just to the acting of U on ip(0) — e/, 



^(0) 



^(1) = Ety(0) 



— tml^m ~t~ t n ie n . 



Since this happens in the same way at all scatterers we 
generalize relation (^|) to the time evolution of an arbi- 
trary state and define 



(9) 



for integer n. With the latter definition the original static 
network becomes a quantum-dynamical model of a dis- 
ordered system. We can relate the eigenvectors 
with eigenvalues {e~ ZUJ } of U to eigenf unctions and en- 
ergy spectrum of a Hamiltonian H by identifying [jl9|2"o 



U = exp(—iHr). 



(10) 



To obtain the LDOS pi(u) = 8(w - u')\<t>f | 2 
(the sum runs over all eigenvalues u) 1 of H) at link I, 
we calculate numerically the temporal Green's function 
G(l,l;t) = (ei\U t ei). In terms of the eigenstates of U the 
Green's function is given by 

G(U;i) = ^|( ei |<m 2 e- Mt = [dup^e-™*. (11) 



Hence the LDOS pi(u>) can be obtained by an inverse 
Fourier transform of G(l,l;t), which on the other hand 
can be generated iteratively according to (0), 
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FIG. 3. Local density of states (LDOS) pi(u>) obtained from 
eq. ([l2|) for a quantum Hall system of 200 x 200 nodes. Due to 
the finite time T (16384 iterations) of the Fourier transform 
each 5-function in the LDOS is broadened into a peak of width 
Alj = 2n/T. 

If 1 f°° 

pdu) = — / dtG(l,l;ty wt = -Re / dtG(l,l;t)e lut . 
2tt J Jo 

(12) 

Fig. (|J) shows as an example the energy dependence of 
the LDOS obtained in that way. 

The three-dimensional network investigated here is 
very similar to the one studied by Chalker and Dohmcn 
. It is built out of layered two-dimensional networks 
with additional inter-layer couplings, that we choose 
slightly different from Chalker and Dohmen pc| . 

In conclusion, we studied numerically the spatial and 
spectral multifractal measures defined by the local den- 
sity of states (LDOS) at the mobility edge of two- and 
three-dimensional disordered electron systems. We have 
presented evidence that both of these measures are equiv- 
alent, as the ratio of their respective generalized dimen- 
sions D(q) and D(q) is simply given by the dimension 
of the system. This result is interpreted as a conse- 
quence of the occurrence of a single energy-dependent 
length scale L u — {p{E c )uj)^ 1 / d acting as the effective 
system size. The numerical calculations were performed 
for two- and three-dimensional network models. These 
model, endowed with a discrete time evolution, turned 
out to be especially suitable for determining the LDOS 
as a function of energy. 
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